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We investigate finite-size effects on the phase structure of chirai and difermion condensates at fi- 
nite temperature and density in the framework of the two-dimensional large- TV Nambu-Jona-Lasinio 
model. We take into account size-dependent effects by making use of zeta-function and compactifi- 
cation methods. The thermodynamic potential and the gap equations for the chirai and difermion 
condensed phases are then derived in the mean-field approximation. Size-dependent critical lines 
separating the different phases are obtained considering anti-periodic boundary conditions for the 
spatial coordinate. 
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I. INTRODUCTION 

Investigation on the phase structure of the strongly 
interacting matter is one of the most interesting top- 
ics in the realm of the standard model for the funda- 
mental forces, in particular for the confining-deconfining 
phase transition. This transition is expected to take 
place far from the asymptotic freedom (high energy) do- 
main of Quantum Chromodynamics (QCD) and so, non- 
perturbative methods are needed in order to approach it 
and other phenomena occurring in this region. The most 
currently used non-perturbative methods involve com- 
puter simulations in lattice field theory and have given 
many interesting results, as for instance numerical esti- 
mates for the confining-deconfining transition tempera- 
ture. In what concerns precise analytical studies, they 
are very difficult in the low energy domain of QCD, due 
to its complex field-theoretical structure. 

In reason of these difficulties, phenomenological mod- 
els for QCD have been adopted along the years. Sim- 
plified effective theories have been largely employed as a 
laboratory to get, analytically, insights on the behavior 
of hadronic matter, particularly in their low-dimensional 
versions. One of the most frequently used of these mod- 
els is the Nambu-Jona-Lasinio (NJL) model It is a 
theory of nucleons and mesons, defined in a spacetime 
with an even number of dimensions, constructed from 
directly interacting Dirac fermions with chirai symme- 
try. Phase transition is considered in a manner analo- 
gous to the appearance of Cooper pairs from electrons 
in the BCS theory of superconductivity. This is the ori- 
gin of the expression "color superconductivity" for the 
hadronic phase transition in the context of this model. 
Nowadays it is well known the usefulness of the NJL 
model in the description of the phase diagram of both 
chirai broken phase (quark-antiquark condensation) and 
color superconducting phase (diquark condensation). In 
addition, the NJL model is specially convenient in the 
study of systems under certain conditions, like finite tem- 
perature 0, H| , finite chemical potential, external gauge 
field, among others 0-0]. 

An interesting aspect in the study of the phase transi- 
tions of the NJL model is the relevance of the fluctuations 



due to finite-size effects in the phase diagram. With this 
purpose, different approaches have been used to study 
various aspects of these effects 0413 • Other phenomeno- 
logical approaches have also been adopted. For instance 
in i2Q] a variational procedure is employed to study finite 
density QCD in a model in which the interaction between 
quarks is supposed to be mediated by instantons. This is 
related to the picture of hadrons as assumed in the MIT 
bag model, where nucleons are considered as droplets 
in a chirally symmetric restored phase. These authors 
find that at densities high enough the chirally symmetric 
phase fills space, and color symmetry is broken by the 
formation of a quark-quark condensate. 

In this paper we intend to generalize, by including fi- 
nite size effects, previous results obtained for finite tem- 
perature 0]. We make use of the techniques intro- 
duced in Refs. [ill Gil an d investigate finite-size effects 
on the phase structure of chirai and difermion conden- 
sates at finite temperature and density in the frame- 
work of the two-dimensional large- A NJL model. This is 
done through zeta-function regularization and compact- 
ification methods [2l| [22[ ■ This approach allows to de- 
termine analytically the size-dependence of the effective 
potential and the gap equation. Then, phase diagrams 
at finite temperature and chemical potential, where the 
symmetric and nonsymmetric phases are separated by 
size-dependent critical lines, are obtained. 

Let us remark that our interest in the two-dimensional 
version of the NJL model is an attempt to investigate the 
qualitative aspects relevant to the chirai and difermion 
condensations under the influence of size finiteness of the 
system. In fact many properties of four-fermion interact- 
ing models are similar in lower and higher dimensions, 
and so, we can expect that results obtained in the 2D NJL 
model reflect properties of a more realistic 4D model. 

It could be argued that spontaneous symmetry break- 
ing does not occur in two dimensions, as a consequence 
of the Mermin- Wagner- Coleman theorem. However this 
is not the case in the present situation. We emphasize 
that although the Mermin- Wagner- Coleman theorem de- 
nies the spontaneous breaking of continuous symmetries 
in two dimensions 12311 . it does not apply in the large- A 
limit @, i, 0, H, [T^nH-ii]. Therefore since our interest 
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is in the analysis of the possible breaking of continuous 
symmetries, it is legitimate to study symmetry breaking 
effects in terms of the low-dimensional NJL model in the 
large- N limit. 

We organize the paper as follows. In Section II, we 
calculate the effective potential of the NJL model in 
the mean-field approximation, using the zeta-function 
method. The analysis of spontaneous symmetry break- 
ing induced by the chiral and difermion condensates are 
also done. The size-dependent gap equations is discussed 
in Section III, while the phase diagrams are shown and 
analyzed in Section IV. Finally, Section V presents some 
concluding remarks. 



II. THE MODEL 

Our starting point is the two-dimensional massless ver- 
sion of the extended NJL model, described by the La- 
grangian density 0, d, [I?} , 

+9D (^ (l) 75^ 0) ) , (1) 

where ip and ip are the fermion fields carrying N flavors 
= l,...,N; repeated flavor indices are summed), fi 
is the chemical potential and the 7 matrices are in the 
two-dimensional space representation, with 7 15 = 7°7 1 . 
Notice that the Lagrangian density possess O(N) flavor 
symmetry and discrete chiral symmetry. In the following, 
unless explicitly stated, we use natural units, K = fcs = 
c = 1. 

Choosing a particular representation we have, 
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where, in this representation, C — — 7 . In this case, the 
pairing term reads 



9d 

9D 



We perform the bosonization by introducing the aux- 
iliary fields er, A, associated to the bilinear forms in 
the above Lagrangian density as: g s 'ip^ip^ 1 ' = and 
gd£y5ipj ipg = A. Therefore the modified Lagrangian 
density becomes, 



C = 



$® (ifi-*- /i7°) ^ ~ \^ (e^'tyW) 
+1 (^i W ^ W ) A " " ^|A| 2 . (4) 



We see from Eq. ([J| that the auxiliary field a plays the 
role of a dynamical fermion mass, such that when it has 



a non-vanishing value, the system is in the chiral bro- 
ken phase. The auxiliary field A is associated with the 
difermion condensate. 

Then, integration over ip and ip^ generates the effective 
action, 
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Notice that the trace over the flavor indices in Eq. §5§ 
gives a factor N, which allow us to set, in the large-iV 
limit, gsN = Gs and goN = Gd, with Gs and Go fixed 
at N — > 00. Thus, the effective potential is obtained in 
the mean-field approximation (i. e., a and |A| uniform) 
from Eq. ©, 



Ueff(cT,A) 



Al : 
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where tr means the trace over spinor indices. 

Our aim is to take into account simultaneously finite- 
temperature and finite-size effects on the phase structure 
of the model; in order to do this, we consider an Eu- 
clidean space, with imaginary time and the spatial co- 
ordinate being compactificd. We denote the Euclidian 
coordinate vectors by xe = (xo,x), with xq £ [0,/3] and 
x £ [0, L], where (3 is the inverse temperature, (i = T^ 1 , 
and L is the size of the system. This corresponds to the 
generalized Matsubara prescription, 
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where 



ko -> w no = ^ ^no + ^ ; no = 0, ±1, ±2, • • • , 
2tt , 

k -> uj n = — (n + c); n = 0, ±1, ±2, 

in the above equation the quantity c is such that, c = 
and c = h for, respectively, periodic and antiperiodic 
spatial boundary conditions. In this article we will re- 
strict ourselves to antiperiodic spatial boundary condi- 
tions, which is a natural choice for fermionic systems. 
The case of periodic spatial boundary conditions would 
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follow along parallel lines. Unless explicitly stated, in 
all cases studied the spatial boundary conditions are an- 
tiperiodic. 

Using Eq. ([S]) we get, after some manipulations, the 
effective potential carrying finite-temperature and finite- 
size effects, omitting terms independent of |A| and <r, 



= ^ 
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where 
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|A| 2 + a 2 + M 2 + W 2 
±2[|A|V + A1 2 K + a 2 )]" 



(9) 
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The effective potential in Eq. © can be rewritten in 
terms of the Epstein zeta- functions, Y(s), that is, 



0*,A) = 



where 
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s=0 
(11) 
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The analysis of the phase diagram of the temperature 
and boundary-dependent model is performed through the 
solutions of the gap equation containing thermal and 
boundary corrections. However, for completeness and 
to set up the free space parameters, in the next section 
we start by treating the zero-temperature model without 
chemical potential and in absence of spatial boundaries. 



III. MODEL AT T = 0, p = AND WITHOUT 
SPATIAL BOUNDARIES 



Let us study the model introduced in the previous sec- 
tion, without compactification of the spatial dimension 
and at zero temperature. This case has been well-studied 
in Ref. Q ; here we perform this study to define the zero- 
temperature free space parameters in absence of chem- 
ical potential. The renormalization conditions for the 
coupling constants are, 
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and 



where Gsr and Gdr are the renormalized coupling con- 
stants, and <7o and A are scale parameters. Then, taking 
/i = T = in Eq. ([8]), and performing the integration 
over pq, the unrenormalizcd effective potential can be 
rewriten as 

a 2 |A| 2 
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where k± is given by Eq. (|10p . with the replacement 
uj n — > k. So,we obtain from (|13[) and (fT4|) the following 
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for £ -4 0. 

Hence now we are able to write the renormalized effec- 
tive potential, 

U eff (a,A) = ai a 2 + a 2 \A\ 2 




(18) 



where 
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(19) 



in the above equation, FP {...} means the finite part of 
the terms between brackets. When we set a.\ = 0, we 
have the vacuum values, a ^ and A = 0. Choosing 
renormalization scales as (Jq — m (m is a scale parameter) 
and Ao = 0, then a 2 > 0, with Gsr. = vr. This is the 
chiral condensate sector. On the other hand, if we take 
the vacuum with A^0 and a = 0, and choosing a 2 = 0, 
ai > 0, Ao = m and <tq = 0, we have Gdr = 2-zr, 
corresponding to the difermion condensate sector. 



IV. MODEL AT FINITE T, (j, AND WITH 
SPATIAL BOUNDARIES 
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Now we take into account temperature, chemical po- 
tential and finite-size dependence. Noticing that the 
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T, fi, L-dependent contributions do not alter the struc- 
ture of ultraviolet divergences discussed in Section IIII1 
then the renormalized effective potential is given by 



otxcr 2 + a 2 |A| 



(20) 



Thus, the T, /z, L-dependent phase diagram of the model 



can be analyzed from U^ff and the gap equations 



A. The chiral condensate sector 

To analyze the pure chiral condensate sector of the 
model, we must take the effective potential in Eq. (|2"0")) 
with A = and a\ = 0, as pointed out at the end of 
Section Uni In order to do this, we perform an analytical 
continuation of the Epstein zeta-function Ha, [HI j 
which gives, 



£o>o*-l] 

+ 2a - cosine) 
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(21) 



Remember that we have fixed m as the scale of the model, 
redefining the relevant quantities as U e f / / m — > U e f / , 
a/m — > a, Lm — > L, fi/m — > /i and /3m — > /3. Firthcr- 
more, we can see that in the bulk limit, L — > oo, Eq. (|21[) 
becomes 
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(22) 



an expression which agrees to that in Ref . Q . 

The ground state is analyzed by means of the minimum 
of the effective potential, which corresponds to the gap 
equation, 
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The non-vanishing solution of the gap equation, yields 
the dynamically generated fermion mass, which comes 



from the equation, 
1 2 °° 

— lncr 2 \^ cos (2nnc)Ko{nLa) 



E E 
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= 0. 



(24) 



B. The difermion condensate sector 

On the other hand, the pure chiral condensate sector 
is studied by taking in Eq. (|2U)) a = and a 2 = 0. Af- 
ter performing the analytical continuation of the Epstein 
zeta-function, Y^, we have the following effective poten- 
tial, 



OA) = ^A>iA*-]] 



2A ^ cos (27rnc) , r . , . . 
> cos (uLn)K 1 (nLA) 
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(25) 



where we have used again m as the scale parameter, 



redefining U e ff/r 



U e ff, A/m -> A, Lm 



/i/m — > /i and /3m — > /3. 

To verify the consistency of the model, we take the 
situation without spatial boundaries, that is, L — > oo. In 
this case, Eq. (|25|) becomes 



m L f {A) = -^A 2 [lnA 2 -l] 
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which is independent of the chemical potential, in agree- 
ment with Ref. [3j]. 

The gap equation for the difermion sector is given by 



(27) 



from which we get the (/3, /i, L)-dependent non-vanishing 
solution in the form, 
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Hence, since A is the order parameter of the difermion 
condensation phase transition, the solution of the gap 
equation (f28| is useful in the characterization of the 
(T, /i, independent phase diagram of the difermion sec- 
tor of the model. 
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FIG. 3: Effective potential for the chiral condensate sector for 
x = 1/L = 0.4 and T = 0.4. Solid, dashed and dotted lines 
represent y, = 0, 0.4 and 0.8, respectively. 



V. PHASE STRUCTURE 



FIG. 1: Effective potential for the chiral condensate sector for 
x — 1/L — 0.2 and T = 0.1. Solid, Dashed and dotted lines 
represent fj, = 0.632, 0.639 and 0.645, respectively. 
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FIG. 2: Effective potential for the chiral condensate sector for 
x = 1/L = 0.3 and T = 0.2. Solid, dashed and dotted lines 
represent fj, — 0.63, 0.64 and 0.65, respectively. 



Now we are able to analyze the L-dependcnt critical 
curves in the phase diagram of the model. First, we 
analyze the effective potential for the chiral condensate 
sector. Remarking that Co is the scale of the theory, 
we can redefine the relevant quantities as Aq/cq Ao, 
Lao L, h/ctq — > /i and /3<7q — > ft. 

In Figs. [T][3] are plotted the effective potential in Eq. 
(|21"|) for different values of fx, x = 1/L and T. Some de- 
tailed information on the influence of spatial boundaries 
can be obtained from these figures; from FigfT] we see 
that a first order transition occurs for increasing values 
of fi, for the values of L and T given (in arbitrary units). 
For a smaller value of L and a larger value of T, (see 
Fig. (2), the phase transition takes place in two steps, as 

is increased: it is a first order transition, but such that 
the absolute minimum of the effective potential is dis- 
placed to a smaller value of a, as /i is increased (fi = 0.63 
and fj, = 0.64 in the figure); for a high enough value 
°f A* (a 1 — 0-65 in the figure), the transition becomes a 
second-order one, when the first non- vanishing extremum 
disappears 31]. Finally, we see from Fig. [3] that there is 
a second order phase transition as \i increases, for large 
enough values of L and T. 

To obtain the phase diagram more precisely in the case 
of a second-order transition, we must determine the val- 
ues of L, T and /x at which the dynamical fermion mass 
vanishes. These critical lines are determined by taking 
the fermion mass equal to zero in the gap equation, 



9 T -jp,L, \ 



0; 



(29) 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
I- 1 

FIG. 4: The phase diagram for the chiral phase transition 
in the (/x, T)-plane. As indicated, the curves are drawn for 
x = 1/L = 0.3,0.4 and 0.5, respectively. The chiral conden- 
sate region is below each curve. Dotted parts of the curves 
correspond to a first-order phase transition, while the solid 
ones represent a second-order transition. 



on the other hand, in the case of a first-order transition, 
we have studied the non- vanishing solution of the gap 
equation , togheter with the behavior of the effective 
potential. 

Taking into account simultaneously both first and 
second-order transitions, the critical lines in the (T, //) 
plane are displayed in Fig. [4] for different values of x = 
1/L. In this figure, solid lines correspond to second-order 
and dotted lines correspond to first-order phase transi- 
tions, respectively. We see that as i = 1/L increases, 
smaller values of temperature are necessary to reach the 
symmetric chiral phase in the region of lower values of 
fx; in addition, in the region of greater values of fj, and 
small T, the chiral broken phase spreads out. In what 
concerns the order of the phase transition, we see that as 
x = 1/L is increased, the region of first-order transition 
is suppressed, occurring only a second-order one. 

Now let us analyze the fermion-fermion pairing phase 
structure. We plot in Figs. [5][7jthc effective potential in 
Eq. (p2ij|) for different values of fi,x = 1/L and T. From 
Fig. [5] we see that for small values of a; = 1/L, that this 
phase transition is of first-order and is independent of 
H in these cases. Fig. [6] shows that for higher value of 
x and smaller temperature, the system undergoes a first 
order phase transition for increasing values of /i. From 
Fig. it can be seen that for fixed /i and a fixed smaller 
temperature, a first-order phase transition occurs as /i 
increases. In these cases the difermion condensate phase 
is destroyed by experencing a first-order phase transition 
as the size decreases or the chemical potential increases. 



0.0 0.2 0.4 0.6 0.8 1.0 

A 

FIG. 5: Effective potential for the difermion condensate sec- 
tor, for x = 1/L = 0.1 and T — 0.294. The lines representing 
fi — 0, 0.5 and 1 coincide. 
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FIG. 6: Effective potential in the difermion condensate sector 
for x = 1/L = 0.5 and T = 0.2. Solid, dashed and dotted 
lines represent [i — 0.700, 0.747 and 0.800, respectively. 

In Fig. [8] is plotted the phase diagram of the system 
in the (T, /x)-plane with different values for x = 1/L. 
From the figure we see that as the size of the system de- 
creases, greater values of the chemical potential are nec- 
essary to reach the difermion condensed phase region. On 
the other hand, for larger sizes (for instance L = (0.1) -1 
in the figure), this region spreads out, and can be reached 
for practically any value of the chemical potential. 



7 



0.211 



0.15 



0.10 



0.05 



0.00 



-0.05 




FIG. 7: Effective potential for the difermion condensate sec- 
tor in the antiperiodic case Solid, dashed and dotted lines 
represent x — 1/L = 0.1,0.49 and 0.6, respectively. 
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FIG. 8: The phase diagram for the difermion condensation 
phase transition in the (fi, T)-plane. The lines represent 
x = 1/L — 0.1,0.35 and 0.5, respectively. The difermion 
condensate region is below each line, which represents a first- 
order phase transition. 



VI. CONCLUDING REMARKS 

In this article we have analysed finite-size effects for a 
version of the NJL model presented originally in Ref. 
for finite temperature and density in free space. This 
has been done by using the extended Matsubara pre- 
scription and zeia-function regularization methods. Fi- 
nite chemical potential effects have also been included. 
These techniques allowed us to construct a renormalized 
effective potential taking into account simultaneously the 
dependence on temperature, the size of the system and 
chemical potential. From this effective potential we have 
obtained the size-dependent gap equations at finite tem- 
perature and density. A throughout analysis of the size 
dependence of the phase diagram has been performed for 
both chiral and difermion sectors. 

In what the chiral condensate sector is concerned, a 
first-order transition occurs for increasing values of the 
chemical potential at fixed values of temperature and of 
the size of the system. For smaller values of L and simul- 
taneously larger temperatures, the system undergoes a 
phase transition in two steps, for increasing values of the 
chemical potential: a first-order transition but such that 
the absolute minimum is displaced to a smaller value of a, 
as /i increases. For some larger value of fi the transition 
becomes a second-order one. We have taken into account 
simultaneously first- and second-order transitions. This 
has been done by a study in the T — fj, plane of the be- 
havior of the system for different values of x — 1/L. We 
conclude that as x increases the region of a first-order 
transition disappears, leaving place for a second-order 
transition. 

The fcrmion-fermion pairing phase structure has also 
been investigated, from an analysis of the effective po- 
tential for different values of fj,, x and T. We conclude 
that for small values of x we have a first-order transi- 
tion, independently of the value of fi. For fixed \x and 
a smaller temperature, a first-order transition occurs as 
x increases. In general, as the size of the system be- 
comes smaller, larger values of fj, are needed to attain the 
difermion condensate sector. 

It should also be noticed that Fig |8] indicates that as 
L decreases there exists a minimal size of the system for 
which the difermion condensate exists at zero chemical 
potential. For smaller values of L, the difermion phase 
reappears if we increase the chemical potential. The full 
phase structure, considering both sectors simultaneously, 
could note be obtained analytically. 
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